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In the study of laser-driven electron acceleration, it has become custoinary to work within the framework of paraxial wave 
optics. Using an exact solution to the Helmholtz equation as well as its paraxial counterpart, we perform numerical simulations 
of electron acceleration with a high-power TMqi beam. For beam waist sizes at which the paraxial approximation was previously 
recognized valid, we highlight significant differences in the angular divergence and energy distribution of the electron bunches 
produced by the exact and the paraxial solutions. Our results demonstrate that extra care has to be taken when working under the 
paraxial approximation in the context of electron acceleration with radially polarized laser beams. © 2013 Optical Society of 
America 
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Electron acceleration driven by high-power lasers has at- 
tracted much effort, since it could lead to the development of 
compact electron accelerators. Such devices are expected to 
be useful in many applications, ranging from coherent X-ray 
generation iHill to electron diffraction imaging ||3]. In laser- 
driven particle acceleration, the challenge is to find an effi- 
cient way to apply the electromagnetic field. Among the many 
proposed methods, the use of radially polarized laser beams 
(RPLBs) for electron acceleration in vacuum has been the 
object of many investigations iH UTTI . In this scheme, which 
was recently demonstrated experimentally lfT2l . electrons ex- 
perience sub-cycle acceleration from the longitudinal electric 
field component at the center of an ultra-intense TMqi beam. 
This mechanism induces a strong longitudinal compression 
effect that could theoretically lead to the production of at- 
tosecond electron bunches ll6l lT3l . 

Many studies that have investigated on-axis EUSJEIID and 
off-axis ||6l[T0l electron acceleration in RPLBs are based on 
the usual paraxial approximation. In these analyses, it is gen- 
erally assumed that the paraxial approximation is valid when- 
ever A:oZr > 125 (or Wo //lo > 2.5), where zr and Wo are respec- 
tively the Rayleigh range and beam waist size, and Aq -Inj ko 
is the dominant wavelength of the laser. 

In this Letter, we question the validity of the paraxial ap- 
proximation in the above-mentioned parameter regime. Using 
a rigorous solution to Helmholtz equation for a TMqi beam as 
well as its paraxial counterpart, we highlight significant dif- 
ferences in the electron acceleration dynamics predicted by 
the exact and the paraxial fields. The exact fields are used to 
investigate the origin and importance of these differences. 

A nonparaxial TMqi beam propagating along the positive 
z axis with its beam waist at z = is described in complex 
notation (i.e., E - Re{£'e-''<"»'"*'), where coq - cko and (po 
are the beam frequency and phase constant) by the following 
field components in cylindrical coordinates (r, (p, z) 111411151 : 



B^ = Boe-*^°";i(M)sin0. 



(3) 



Er = -jEoe''">"j2{koR) sin cos 



E, = -jjEoe-'"'" [joikoR) + j2(koR)P2{cosd)] 



(1) 
(2) 



Here Eq - cBq are amplitude parameters, ^o is the beam 
wavenumber, R - [r^ + (z + y'fl)^]'^^ is the complex radius, 
COS0 = (z + ja)/R defines the complex angle 0, j„{-) is the 
«-th order spherical Bessel function of the first kind, and 
P2( ) is the Legendre polynomial of degree 2. The parame- 
ter a is a real and positive constant called the confocal pa- 
rameter. The latter may be used to characterize the degree 
of paraxiality of the beam since it is monotonically related 
to the Rayleigh range and beam waist size by the relation 
Zr - kow^/l - [ -sjl + (koa)^ - l]/ko. The power carried along 
the z axis by a nonparaxial TMqi beam is lfT6l 



Poe 



kla^ 



\^kQa sinh(2A;oa) -cosh(2^o'2) + 1 -2A:Qa^j , (4) 



where Po = AEci?' I^Tlak^- Note that the fields described in 
Eqs. ([T]i-(l3]i represent a rigorous solution to Helmholtz equa- 
tion. We will thus refer to them as the exact TMoi fields. 

In the paraxial limit, namely when k^a » 1, the fields ([T]i- 
^ can be expanded as power series of the parameter 6 - 
l/ki)a. Using the normalized coordinates p - r(2a/A:o)"'^^ 
and ( - zja (since zr - kow^/l » a), we find, up to terms of 
order S^: 



E, = -^_Eo[pf-6y' - {3jpf'+3p'f'-ip'f')6''' 

E: = jEo[if-jp'f)S' 

- ijf+5pY - '4p'f - kp'f)S' + 0(6")] oxpij^') 
B, = -^Bo[p/^^^^ - Upf'+2p'f'-ip'f')6'l' 



(5) 



(6) 



(7) 



where / = l/(f + j) and - -{(/6 + p^f). The lowest-order 
terms of Eqs. (I5]l-(|7]i correspond to the well-known paraxial 
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TMol fields, which are commonly used to analyze electron 
acceleration in RPLBs I 



f(0) 



Ef^ ^jEoif-jp-f)6'-cxp{j^') 



(8) 
(9) 



To simulate electron acceleration in RPLBs, we numeri- 
cally integrate the Newton-Lorentz equations: 



dr 

dt 



d\ 
dt 



ynie 



E + vxB- — (vE) 



(10) 



where e, nie, r, v are the electron's charge, mass, position, ve- 
locity, respectively, and y = (1 - |vp/c^)"'^^. We also suppose 
that the laser beam is pulsed, which we model by multiply- 
ing E and B by sech(^/^o)^ where ^ = wof - fcoZ- This is to 
ensure that the fields satisfy Maxwell's equation in the limit 
^0 » 1 for any value of the phase ^ 121. Furthermore, we use 
Aq - 800 nm, although it can be readily shown that the results 
are scalable to any value of Aq Igl fTTI . 

We consider a cloud of electrons initially at rest in the 
(r, z) plane outside the laser pulse. The initial position of the 
electrons are drawn randomly from a two-dimensional Gaus- 
sian distribution centered at the origin with standard deviation 
cr,. = cr. - /lo/lO. The electrons are accelerated by a pulsed 
TMoi beam with k^a - 500 (which corresponds to wo ~ 5Aq), 
a value generally considered well inside the paraxial regime. 
Space-charge eff'ects are neglected; each trajectory is com- 
puted independently. The relativistic electron bunch produced 
after performing the numerical simulation with the exact and 
the paraxial TMoi fields from the same initial conditions are 
shown in Fig. [T] A comparison of the datasets in Fig. [11 a)- 
(c) immediately shows an enormous difference in the bunch 
transverse extent; the angular divergence of the bunch ac- 
celerated by the exact fields is about 100 times larger than 
the bunch accelerated by the paraxial fields. Moreover, the 
electron bunch energy gain distribution is very different from 
one case to the other, as shown in Fig. [lid). The average en- 
ergy gain obtained with the paraxial fields is however close to 
the energy gain near the optical axis with the exact fields, as 
shown in Fig. [He). This agrees with previous results reported 
in lim for on-axis acceleration. Since the near-axis electrody- 
namics is similar in both cases, the strong longitudinal com- 
pression predicted by the paraxial fields is also observed with 
the exact fields. 

To understand why such important discrepancies between 
the exact and the paraxial cases arise for ofF-axis electrons, it 
is instructive to look at the equations of motion. When elec- 
trons close to the optical axis interact with the laser beam, 
they are primarily accelerated in the positive z direction by 
the longitudinal electric field. For an electron with <s: v,, 
the equation governing the radial velocity is approximately 



dVr 
It 



(11) 



From Eq. ((8]l, we see that the paraxial E^ and field com- 



ponents are perfectly symmetric, i.e., = cBj^"'. Therefore, 
according to Eq. (fTTT i. an electron travelling primarily along 
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Fig. 1 . (Color online) Electron bunch accelerated by a TMqi 
pulsed beam with P = 10'^ W, k^a = 500, = 14.21, 
00 = ^- The simulation was performed with: (i) the exact 
fields [Eqs. Q-Q], (ii) the paraxial fields [Eqs. ([Ill-®], and 
(iii) the corrected paraxial fields [Eqs. ([5}-([2) up to 0{5^)\. (a) 
Bunch snapshot and (b)-(c) close-up view of the front end; 
(d) energy gain distribution; (e) average energy gain versus fi- 
nal radial coordinate. The electrons are initially at rest outside 
the laser pulse in the (r, z) plane and distributed randomly ac- 
cording to a Gaussian distribution centered at the origin with 
- CTj. = AqIIQ. The results are computed 15 ps after the 
passage of the beam at z = 0. Only = 200 electrons are 
shown in (a)-(c), while = 50000 different initial conditions 
were used to obtain the results in (d) and (e). 



the optical axis at a relativistic velocity will feel a quasi-null 
force in the radial direction. This explains the observation 
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Fig. 2. (Color online) Normalized amplitude of (a) —eE- and 
(b) -e(Er - cBij,) at f = (exact fields). During sub-cycle ac- 
celeration, electrons will spend most of their time between a 
maximum of -eE. and the minimum located behind [an ex- 
ample is shown by the rectangle in (a) and (b)]. In this phase 
configuration, electrons travel at V; a; c and feel, according to 
(b), a radial force directed outward. Beam parameters are the 
same as in Fig.[I] 

of electron bunches with very narrow transverse extent when 
performing the simulations with the paraxial fields. 

However, this argument does not hold in the exact case. Ac- 
cording to Eqs. ([TJ and Q, the exact Er and field compo- 
nents are not perfectly symmetric. The perturbative series (|5]i 
and (|7]i show that the symmetry between Er and is broken 
as soon as the first nonparaxial correction is introduced. Since 
Er and B^ reach their maximum value at p » [(1 + ^^)/2]'''^, 
the ratio of the maximum value of the radial field term in 
Eq. (fTTT i. \Er - cB^I, and the accelerating field, l^-l, scales as 
5*^^(1 + f^)'^^. The radial force component may therefore not 
be neglected when koa > 100, and its relative importance be- 
comes greater as ( increases. Moreover, Fig. |2] demonstrates 
that the phase configuration in which electrons spend most of 
their time during sub-cycle acceleration is associated with a 
radial force directed outward. This explains why the angular 
divergence of the bunch accelerated by the exact solution is 
larger instead of smaller 

To verify explicitly if the discrepancy between the electron 
bunches in Fig.[T]originates from the articifial symmetry be- 
tween the paraxial Ef' and field components, we have 
performed the same numerical simulation with the corrected 
paraxial fields, i.e., the paraxial fields with the first nonparax- 
ial correction to each component [Eqs. (|5)-(l7) up to 0{6^)]. 
The results obtained are shown in Fig.[T] We immediately see 
a much better agreement with the exact case. While the final 
positions of the electrons do not match perfectly, the spatial 
configuration of the bunches and their energy gain distribu- 
tion are much more similar Note that the small diff'erences in 
the energy gain distribution can be understood from Fig. [lie). 
Indeed, we see that the corrected paraxial results get less ac- 
curate as p increases. This is a consequence of the fact that 
the magnitude of the nonparaxial corrections increases as we 
move further from the optical axis. Taking into account the 
first nonparaxial correction to Er and B^ thus only offers a 
limited solution for electrons far from the optical axis. 

Despite the fact that the validity of the paraxial approx- 
imation was never fully adressed in the literature for k^ZR 
values well above 100, an order of magnitude comparison 



may be performed with existing results. In ifTSl . an RPLB 
with Wo = 3/lo and second-order field accuracy was used to 
accelerate electrons from a target of size comparable to the 
initial conditions used here. The resulting electron bunch an- 
gular divergence was reported to be approximately A0 ^ 3°. 
With the parameters used in this Letter, we obtain AQ ^ 5° 
with the exact fields, which is of the same order of magnitude 
as in ifTJl . compared to < 0.1° for the paraxial fields. 

In conclusion, in a parameter regime where the paraxial ap- 
proximation was previously considered valid, we have high- 
lighted significant differences between the properties of elec- 
tron bunches accelerated by paraxial and exact TMoi beams. 
These differences originate from the symmetry between the 
paraxial Ef^ and fields. This artificial symmetry is bro- 
ken as soon as the first nonparaxial corrections to the electro- 
magnetic field are taken into account, which allows to obtain 
more accurate results. The considerations presented in this 
Letter are also believed to apply to ultrashort pulses, since 
the relation - cB^ always holds for the paraxial TMqi 
fields, regardless of the pulse duration Js]. Our study thus ad- 
vocates that special care has to be taken when working under 
the paraxial approximation in the context of electron accel- 
eration in RPLBs. It should be reminded that under relativis- 
tic conditions, nonparaxial field corrections may always yield 
major differences in the trajectories of ofF-axis electrons, even 
for very large values of k^a. Solutions as accurate as possible, 
ideally exact, should thus be used. 
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